Feshbach resonances with large background scattering length: interplay with 

open-channel resonances 
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Feshbach resonances are commonly described by a single-resonance Feshbach model, and open- 
channel resonances are not taken into account explicitly. However, an open-channel resonance near 
threshold limits the range of validity of this model. Such a situation exists when the background 
scattering length is much larger than the range of the interatomic potential. The open-channel 
resonance introduces strong threshold effects not included in the single-resonance description. We 
derive an easy-to-use analytical model that takes into account both the Feshbach resonance and 
the open-channel resonance. We apply our model to 85 Rb, which has a large background scattering 
length, and show that the agreement with coupled-channels calculations is excellent. The model can 
be readily applied to other atomic systems with a large background scattering length, such as 6 Li 
and 133 Cs. Our approach provides full insight into the underlying physics of the interplay between 
open-channel (or potential) resonances and Feshbach resonances. 

PACS numbers: 34.50.-s, 02.30.Mv, 03.75.Nt, 11.55.Bq 
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I. INTRODUCTION 

A magnetically induced Feshbach resonance is an 
indispensable tool to control the atom-atom interaction 
in ultracold gases. By simply changing the magnetic field 
around resonance, the s-wave scattering length, which is 
a measure for the strength of the interactions, can be 
given basically any value. Exactly on resonance, the scat- 
tering length is infinite, and its value is therefore much 
larger than any other lengthscale that characterizes the 
atomic gas system. This means that the scattering length 
effectively drops out of the physical problem. For very 
low temperatures and on resonance, universal behavior 
has been predicted 0, For higher temperatures, or 
for situations further away from resonance, knowledge 
about the correct energy-dependence of the scattering 
phase shift will be needed to account for an accurate de- 
scription of a resonant two-body interaction |5j, |(j . 

Crossing the Feshbach resonance, the scattering length 
a goes from positive to negative values through infinity, 
and the effective interaction changes accordingly from 
repulsive (a > 0) to attractive (a < 0). Using two atomic 
spin states of a fermionic gas, it is possible to study the 
crossover region between a BCS superfluid and a Bose- 
Einstein condensate (BEC) of molecules @,H,E|- 

The Feshbach resonance is associated with a molecu- 
lar state just below the collision threshold. When the 
magnetic field is changed during an exp eriment, (quasi- 
bound) molecules can be formed LflUJ, 121 while sweep- 
ing through the resonance. At JILA, coherent oscillations 
between atoms and molecules have been observed by ap- 
plying two mag netic field 'Ramsey' pulses close to res- 
onance [l3L The oscillation frequency was directly 
related to the binding energy of the molecules. More 
recently, Feshbach resonances were utilized in the for- 
mation of ultracold molecules from a degenerate atomic 
Fermi gas 0, 0, 0, 0] . Here molecules could be pro- 



duced reversibly by sweeping twice through resonance. In 
some cases these dimers have been cooled further down 
below the BEC transition temperature 0, 0] • In 
the case of bosonic systems, Feshbach resonances have 
been utilized to produce ultracold sodium, cesium and 
rubidium molecules [51 12I l24| . 

Close to the magnetic field value of resonance Bq, the s- 
wave scattering length a shows a characteristic dispersive 
behavior and is given by 



a(B) = a bg 1 - 



AB 

B-Br 



(1) 



The background part of the scattering length ab g sum- 
marizes the effect of the direct scattering processes in the 
open channel, without coupling to other closed spin chan- 
nels. AB is the field width of the Feshbach resonance. 
An example of this behavior can be seen in Fig. ^ Since 
the Feshbach resonance is responsible for the dramatic 
behavior of the scattering length, the background pro- 
cess is usually considered not very interesting, and is in 
some cases even neglected. The direct process is then 
associated with non-resonant scattering, where the value 
of abg should be of order of the van der Waals potential 
range r = (yuCg/Sft 2 ) 1 / 4 [2j|, with fi the reduced mass 
and Cq the van der Waals coefficient. This, however, is 
only true when there is no potential resonance close to 
the collision threshold. 

The effect of potential resonances can be easily esti- 
mated by comparing the value of ctb g to the potential 
range rrj. For example, in the case of 87 Rb, the direct in- 
teraction potential is non-resonant since ab g = 100 ao |26| 
and rrj is of the same order for rubidium. There are also 
situations where this is not the case. When ab g is large 
and positive, there is an open-channel bound state very 
close to threshold. For ex amp le, this occurs for 133 Cs, 
where db g = 905 ao [13, l28| is much larger than r . 
In other situations ab g can be large and negative, in- 
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FIG. 1: Scattering length as a function of the magnetic field 
for 87 Rb in the \f,mf) = |1, 1) hyperfine channel. The 
horizontal dashed line indicates the direct scattering length 
ab g = 100 ao (with ao the Bohr radius) , and the vertical dot- 
ted line indicates the resonant magnetic field Bo = 1007.4 G. 
The field width is given by AB = 0.2 G. 



dicating that the interaction potential has a virtual state 
near threshold. Examples here are 85 Rb with = 
-443 a [3, and 6 Li where a hg w -2000 a [M Ell- 
in these cases, the non-resonant description of the back- 
ground scattering process is not correct. This can be seen 
from a comparison between the single-resonance Fesh- 
bach model and a full coupled-channels calculation over 
the range of energies of interest for current experiments, 
as will be shown in the following. 

In this paper, we demonstrate how the potential reso- 
nance can be taken into account analytically, by properly 
describing the interplay between the open-channel reso- 
nance and the Feshbach resonance. The open-channel 
resonance is associated with a pole of the direct part of 
the scattering matrix. In our approach, we include this 
open-channel pole into the multi-channel Feshbach for- 
malism. The resonance phenomena discussed here in the 
context of cold atomic collisions are quite general, and 
have been studied before in the context of nuclear physics 
[3l|. electron-atom collisions [HI, and electron-molecule 
collisions [3j| • In previous work the potential resonance is 
indirectly included in the closed-channel subspace, how- 
ever, the potential resonance is actually located in the 
open-channel subspace. In our approach, the distinction 
between open-channel and closed-channel resonances is 
clearly made, which results in a better understanding of 
the underlying physics. 

The paper is outlined as follows. In Sec. [D] we intro- 
duce the basic ideas behind our approach. In Sec. lIIII wc 
briefly discuss the theory of Feshbach resonances, based 
on a projection-operator approach. In Sec. IIVI we dis- 
cuss the properties of potential resonances, and explain 
how the so-called Mittag-Leffler expansion can be used 



to account for these resonances. The Mittag-Leffler ex- 
pansion is used in Sec. \V\ to discuss the interplay be- 
tween Feshbach resonances and potential resonances. In 
Sec. IVII we give a thorough discussion of the results of 
our model. We compare our model with other commonly 
used approaches in Sec. IV 111 and conclusions are drawn 
in Sec. EES 



II. BASIC IDEAS 

In this section, we introduce the basic ideas which are 
needed to describe the interplay between potential res- 
onances and Feshbach resonances. We first introduce 
the Feshbach resonance with a non-resonant background 
scattering process, i.e. with a background scattering 
length on the order of the range of the interatomic poten- 
tial. Then we discuss a single-channel s-wave resonance, 
and its relation to the poles of the scattering matrix. 
In case the background scattering length is large, this 
indicates that the open channel is nearly resonant, even 
without coupling to other channels. Then, we investigate 
the effect of the open-channel resonance on the Feshbach 
resonance. We study two different cases: one with large 
positive and one with large negative a^g. 

Feshbach resonances 

The location of the vibrational bound states in the in- 
teraction potential determines the scattering properties 
in ultracold atomic collisions. In a multi-channel collision 
the incoming channel may be coupled to other (open and 
closed) channels during the collision. A channel is ener- 
getically open (closed) when the total energy of the two- 
atom system is above (below) the dissociation threshold 
energy of the incoming channel. In the Feshbach for- 
malism the total Hilbert space describing the spatial and 
spin degrees of freedom is divided into two subspaces, in- 
dicated by V and Q. Generally, the open channels are 
located in P-space, and the closed channels in Q-space. 
As a model example, we consider the case with only one 
open (P), and one closed (Q) channel. A schematic illus- 
tration of the interaction potentials associated with these 
two channels is shown in the left part of Fig- El Since the 
total energy is above the P-threshold, this P-channel is 
automatically the incoming and outgoing channel. Later 
on, we also consider the case where the total energy is 
below the P-threshold. Strictly speaking, the P-channel 
is then energetically closed and should belong to the Q 
subspace. However, for convenience we will still use the 
labelling P and Q to distinguish between these two chan- 
nels. 

In the case of alkali atoms, the P and Q channels are 
associated with different spin configurations, and have a 
different effective magnetic moment. The energy differ- 
ence between the channels can be tuned by changing the 
external magnetic field. This means that the P and Q po- 
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FIG. 2: (Color online) Left part of figure: Schematic illustra- 
tion of the P-channel (dashed) and Q-channel (solid) interac- 
tion potential for B = 0. Right part of figure: The interac- 
tion potentials are asymptotically connected to the magnetic 
field dependent two-particle hyperfine-Zeeman eigenenergies. 
Also shown are the highest P-channel vibrational bound state 
(fp.max) and a Q-channel vibrational bound state (vq) close 
to the collision threshold. At the resonant magnetic field Bo 
the coupling between the Q-channel bound state and the P- 
channel scattering state becomes resonant. 

tentials can be shifted relative to each other. The Zeeman 
energy shift of the P-channel and Q-channel thresholds is 
schematically shown in the right part of the figure. Also 
indicated is the energy of the highest P-channel vibra- 
tional bound state (ypmax) an d a Q-channel vibrational 
bound state {vq) close to the P-threshold. 

During the collision the atoms are coupled from the 
P-channel to the Q-channel. The scattering process be- 
comes resonant when a Q-channel bound state is located 
close to the P-channel collision threshold, giving rise to a 
Feshbach resonance. The unperturbed Q-channel bound 
state is dressed by the coupling to the P-channel. This 
dressed state can be considered as a (quasi-)bound state 
of the total scattering system. The colliding atoms are 
temporarily captured in the (quasi-)bound state, and 
after a characteristic time r = 2h/T, with T the de- 
cay width of the (quasi-)bound state, they return to 
the P-channel. At the magnetic field value Bo, where 
the dressed state crosses the P-threshold, the scattering 
length has a singularity. The resulting behavior of the 
scattering length was already shown in Fig. ^ 

The scattering length is related to the zero-energy limit 
of the scattering matrix. In the Feshbach description the 
energy-dependent scattering matrix is divided into two 
parts. The direct part S^ rcct describes the scattering 
processes in the P-channel, without coupling to the Q- 
channel. The Feshbach resonance is described by the 
resonant part 5^ s . The resulting total scattering matrix 
S = Sd ilect S® s summarizes the energy and field depen- 
dence of the scattering process around resonance. 

The direct part of the scattering process is usually 
assumed to be non-resonant. In the treatment by Mo- 
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FIG. 3: Schematic illustration of the poles of the S matrix 
in the complex k plane. The bound-state pole is located at 
the positive imaginary k axis. For a less attractive potential, 
this pole will move down the imaginary fc-axis, giving rise to 
a virtual state. 

erdijk et al. [34[ the direct part of the scattering ma- 
trix is described by a single parameter only: S^ ilect (k) = 
exp[— 2i/cabg], where k is the relative wavenumber of the 
two colliding particles. In the case Ob g is of order of the 
range of the potential ro, the resulting single-resonance 
Feshbach model accurately summarizes the scattering 
process for the energies and magnetic fields of interest. 



Potential resonances 

The direct scattering properties are determined by the 
location of the bound states in the P-channel potential. 
In case the P-channel has a bound state close to thresh- 
old, this channel is nearly resonant and ab g is large and 
positive. This has several important implications, which 
have not been considered before in the description of 
Feshbach resonances in cold atomic collisions. The P- 
channel potential resonance introduces a resonance en- 
ergy dependence of S^ iTect , not included in the simple 
description based on db g only. It is not straightforward 
to describe this resonance by means of a projection for- 
malism analogous to the Feshbach approach, since the 
bound state is located in the open channel. Still, it is 
possible to understand the resonance behavior by study- 
ing the poles of the scattering matrix. 

There is a general connection between the resonances 
and bound states of a scattering system, and the various 
poles of the scattering matrix [33. For single-channel s- 
wave collisions, the poles in the upper half of the complex 
fc-plane can only be located on the imaginary fc-axis. Suf- 
ficiently close to the origin (k = 0), the poles in the lower 
halfplane can only be located on the imaginary axis as 
well [36|. As a result, the energy E — h 2 k 2 /2/i associated 
with these poles will be strictly negative. The P-channel 
bound state is related to a pole of S^ iTect on the positive 
imaginary /c-axis, as indicated in Fig. [3J If the bound 
state is located just below threshold, this pole is located 
close to the origin. 
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Now imagine that the potential gradually becomes less 
attractive. Then, the bound states will move towards 
threshold, and the poles of the S matrix will move down 
the imaginary fc-axis. Consequently, the bound-state 
pole close to the origin will move towards k — 0, cross the 
origin when the bound state is located exactly at thresh- 
old, and eventually turn into a so-called virtual-state pole 
on the negative imaginary fc-axis. The energy associated 
with this virtual state is negative, but there is no proper 
physical bound state associated with this energy. A vir- 
tual state can be regarded as a nearly-bound state that 
behaves much like a real bound state in the inner region 
of the interaction potential. Only in the asymptotic re- 
gion (r — ► oo ) the virtual state 'discovers' it does not 
quite fit to the size of the interaction potential, and the 
virtual state exponentially explodes. 

The direct part of the scattering matrix can be di- 
vided into a resonant and a non-resonant part: S^ iTect = 
S^ g S^ cs . The background part is now related to the true 
range of the potential, while the resonant part takes the 
form: 



k + in 
k — ik 



(2) 



where in is the location of the pole on the imaginary 
/c-axis. In the zero-energy limit, the direct part of the 
scattering matrix is related to the background scatter- 



ing length as ab g = a£ + 1/k. 



Here is of order 



ro, while the resonant contribution of the P-channel is 
given by 1/k. The important effect of a bound or virtual 
state near threshold is the introduction of a resonance 
energy dependence of <5^ rect . More details will be given 
inSecHVl 



Interplay 

Now we return to a multi-channel collision, with a 
Feshbach resonance resulting from coupling to Q-space, 
and a potential resonance associated with P-space. The 
total scattering matrix takes the form S = S^ g S^, s S^, s , 
where the background part is related to the true range of 
the potential. In a multi-channel problem, the bound and 
quasi-bound states of the scattering system are related to 
the poles of the total S matrix. The bound-state poles 
are located on the positive imaginary fc-axis, whereas the 
quasi-bound state poles are located in the lower half of 
the complex fc-plane. These poles are generally not lo- 
cated on the imaginary axis, as they are related to com- 
plex energies E^ — iT/2. Here Pr is the energy, and T the 
decay width of the quasi-bound state. The decay width 
is related to the lifetime of the quasi-bound Feshbach 
resonance state. 

In Fig.^J we consider the case where the P-channel has 
a bound state just below threshold. This means that ab g 
will be large and positive, which is the case for 133 Cs, for 
example. The bare bound states of the P and Q channels 
are indicated by the horizontal (i/p !inM ) and slanted (vq) 




Magnetic field (arb. units) 

FIG. 4: (Color online) Example of an avoided crossing be- 
tween a P-channel and a Q-channel bound state. The energy 
of the bare bound states is indicated by the horizontal and 
slanted lines. The energies of the dressed states are indicated 
by the solid lines. One of the dressed states crosses the P- 
threshold at Bo- The bare Q-channel bound state crosses the 
P-threshold at Bo- Note that the energies are give relative to 
the P-channel collision threshold (horizontal solid line). 



lines, respectively. These bare states are eigenstates in 
the uncoupled V and Q subspaces. The bare Q-channel 
bound state crosses the P-channel threshold at the mag- 
netic field value Bq. The dressed states are indicated by 
the solid lines. 

Below the P-threshold, the behavior of the dressed 
states resembles an avoided crossing as described by a 
two-level Landau-Zener model. In the inset, the behav- 
ior of the dressed state is shown close to the P-channel 
threshold. At B one of the dressed bound states crosses 
the collision threshold, acquires a decay width, and turns 
into a quasi-bound state. Exactly at Bq, the s-wave 
scattering length goes through infinity. Note that the 
binding energy of the dressed state curves quadratically 
towards threshold according to the well-known relation 
e bind (P) = -h 2 /(2[ia 2 (B)) oc (B - P ) 2 , as follows from 
Eq. . At threshold, the slope of the dressed state en- 
ergy shows a discontinuity resulting from a non-zero de- 
cay width above threshold. We stress that, although the 
behavior of the dressed states shows some resemblance 
with the two-level Landau-Zener description, this model 
does not include these threshold effects and can not be 
used to properly describe the interplay between a poten- 
tial resonance and a Feshbach resonance. 

For dbg large and negative, there is a P-channel vir- 
tual state near threshold. Examples here are 85 Rb and 
6 Li. An example of the energy of the dressed states for 
this situation is given in Fig. [SJ Here the difference with 
the Landau-Zener model is even more striking. The be- 



5 




Magnetic field (arb. units) 



FIG. 5: (Color online) Effect of virtual state on Feshbach 
resonance in 85 Rb in the |2, —2) hyperfine channel. Shown 
are the energy of the unperturbed Q-channel bound state 
(dotted line), and of the (quasi-)bound state (solid line) 
which is 'dressed' by the coupling to the P-channel. The 
(quasi-)bound state crosses the collision threshold at Bo, 
the unperturbed Q-channel bound state crosses the collision 
threshold at Bq. Note that the energies are give relative to 
the P-channel collision threshold (horizontal solid line). 



havior of the dressed state energies can not be under- 
stood on the basis of an avoided crossing between the Q- 
channel bound state and the P-channel threshold. Based 
on a Landau-Zener approach, the dressed state above 
P-threshold would be located at the righthand side of 
the bare Q-channel bound state, which clearly is not the 
case. The binding energy still curves quadratically to- 
wards threshold, which it crosses at the magnetic field 
Bq. From the inset it can be seen that the slope of the 
dressed state energy shows a discontinuity at threshold, 
which is even more distinct than the situation given in 
Fig. 0] For larger collision energies, the dressed quasi- 
bound state curves back to the bare Q-channel bound 
state [3|j. 

In the following we derive an analytical model that de- 
scribes the interplay between potential resonances and 
Feshbach resonances. We used this model to obtain 
Figs. 01 and |31 which are model examples applicable to 
133 Cs and 85 Rb, respectively. In Sec. lVII we show that the 
agreement between our model and full coupled-channels 
calculations is excellent for a large domain of energies 
and magnetic fields, while the single-resonance Feshbach 
model gives poor results in the same domain. We com- 
pare our model with other approaches in Sec. IVIII 



III. FESHBACH RESONANCE THEORY 

Feshbach resonances in two-body collisions are related 
to the coupling of different spin channels, and can be 
conveniently described in an approach due to Feshbach 
PI ■ In this approach the total Hilbert space H de- 
scribing the spatial and spin degrees of freedom is divided 
into two subspaces V and Q. In general V contains the 
open channels and Q the closed channels. The S and T 
matrices, which are related to the transition probabili- 
ties of the scattering process, are separated in two parts 
accordingly. The P-part describes the direct interactions 
in the open-channel subspace, and the Q-part describes 
the effect of the coupling to the closed-channels. Usually 
the Q-part contains the resonances, and the P-part is as- 
sumed to be non-resonant. We discuss in Sec. lIVI whv in 
some cases the P-part can introduce resonant features as 
well. In Sec. we show how these P-channel resonances 
can be taken into account analytically. First we discuss 
the Feshbach projection formalism in this section. 

One can construct projection operators P and Q, 
which project onto the subspaces V and Q, respectively. 
The Schrodinger equation for the two-body collision can 
then be written as a set of coupled equations: 

(E tot -H PP )\^p) = H PQ \^ Q ), (3) 
(E tot -H QQ )\* Q ) = H QP \^ P ). (4) 

Here we use the notation \^ P ) = P|*), |*q) = 
H PP = PHP, Hpq = PHQ, etc., and E tot = E + E thl 
is the total energy of the colliding atoms, with E the 
kinetic energy and -Ethr the energy of the open-channel 
threshold. As already mentioned, we are only interested 
in scattering processes with only one open channel. The 
P-channel is then simultaneously the incoming and out- 
going channel. Also note that all energies in the follow- 
ing are given with respect to the open-channel collision 
threshold. 

We multiply Eq. (0} from the left with the resolvent 
(or Green's) operator Gqq(E + ) = [E + — Hqq] 



-l. 



I*Q> 



1 



E+ -H, 



QQ 



-Hq P \Vp), 



(5) 



where E + = E + i5 with 6 approaching zero from positive 
values. Substituting the expression for \*&q) in Eq. 
the problem in the V subspace is equivalent to solving 
the Schrodinger equation (E — H c s)\^ P ) = 0, where the 
effective Hamiltonian is given by 



Heft = Hpp + HpQ 



1 



E+ - H QQ 



-H, 



QP- 



(6) 



The first term in the effective Hamiltonian describes the 
direct effect of the open-channel subspace V on the scat- 
tering process. The second term in the effective Hamilto- 
nian describes the coupling of V space to Q space, prop- 
agation through Q space, and coupling back to V space 
again. 
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The resolvent operator Gqq can be expanded in terms 
of the discrete (|<fo)) and continuum (|<^>(e))) eigenstates 
of H QQ : 



Hr 



E 



E-S 



E+ -e 



de. 



(7) 



Now suppose there is only one discrete bound state of 
Hqq for which the eigenvalue e£ is close to the colli- 
sion energy E of the scattering state. That is, the other 
Q-channel bound states are located sufficiently far away 
from the P-threshold so that their contribution can be 
safely neglected. In this case we can also neglect the 
continuum expansion of Hqq, and the problem in the V 
subspace reduces to 



(E-H PP )\^ P ) = H PQ 



E - 



-qHq P \^ P 



(8) 



Now we can formally solve the coupled problem by 
multiplying from the left with the resolvent operator 
G PP (E+) = [E+ -Hpp]' 1 : 



H 



pp 



Hpq^Hq P \* p) . (9) 



The scattering state \^ P ) corresponds to the homoge- 
neous part of the Schrodinger equation projected into 
V space, and is an eigensolution of H PP with outgo- 
ing spherical wave boundary conditions. The scatter- 
ing states are energy-normalized as (^p~(E)\%>p~(E')) — 
S(E — E'), where 5(E) is the Dirac delta function. 

Before calculating the S and T matrices, we introduce 
the Lipmann-Schwinger equations for |^p), where the 
superscript + (-) indicates outgoing (incoming) spherical 
wave boundary conditions [4l|: 



1 



I*?} = \xp) + E± _ Hpp v pp\xp) 



(10) 



Here the Hamiltonian in V space is written as H pp = 
H PP +V PP , where H PP represents the sum of the relative 
kinetic energy operator and the two-particle hypcrfinc- 
Zeeman interactions, and V pp the interatomic interac- 
tions, both projected onto the P-channel subspace. The 
unscattered states \xp) are eigenstates of H pp . 

The Tp matrix, giving the transition amplitude due to 
scattering in the V subspace only, can now be calculated 
as 



(Xp\Tp\xp) = (xp\Vpp 



E+ - Hpp 



Vpp 



(Xp\Vpp\*$). 



\xp) 

(11) 



Introducing the second term in the effective Hamiltonian 
Eq. according to the two-potential theorem |4lJ the 



total T matrix is given by 



(Xp\T\Xp) = (Xp 

(*p\Hpq 



Vpp + Hpq 



E+ - H, 



-H 



QQ 



QP 



I*p) 



<Xp\Vpp\*$) 



E+ -H, 



QQ 



-H QP \*p), (12) 



where we have used the formal solution for \^p) given 
by Eq. © together with relation Eq. (|TU|l . We see that 
apart from the direct term Tp, the transition amplitude 
contains a term which results from the coupling to the 
closed-channel subspace Q. 

The relation between the S matrix S(k) — (xp\S\Xp) 
and the T matrix T(k) = (xp\T\xp) IS given by the oper- 
ator equation 5 = 1 — 2iriT. From this equality, together 
with the expansion of the resolvent Gqq, it follows that 
the S matrix of the effective problem in V space is given 
by 



S(k) = S P (k) -2m 



(^p\HpQ\4> b )(<t> b \HQp\^p) 

E-S 



(13) 



Here S P (k) = (^!p(E)\^+(E)) is the direct part of the S 
matrix, which describes the effect of the scattering pro- 
cess in V space only, without coupling to the closed chan- 
nels in Q space. 

We can finally solve for \^fp) by multiplying Eq. 
from the left with {4>i,\Hqp, which results in 



S(k) = S P (k) 1 - 2m- 



k\H Q p\^] 



E 



A(E) 



The term 



A(E) = (<j> b \H 



QP 



E~^~ — Hp p 



-H 



PQ\ 



(14) 



(15) 



in the denominator is the complex energy shift, which is 
the energy difference between the bare bound state \4> b ) 
and the dressed (quasi-)bound state. 

Inserting a complete set of eigenstates of Hpp, the 
complex energy shift can be written as A(E) = A ICS (E) — 
~T(E). The real part A rcB (E) shifts the unperturbed 
energy e?, and the imaginary part T(E) turns the un- 
perturbed bound state \<p b ) into a quasi-bound state. 
Note that A(E) is purely real for energies below the P- 
threshold (i.e., E < 0). The energy of the dressed states 
can be found by solving for the poles of the 5* matrix. 

In the low-energy domain it is usually assumed that 
the real part of the energy shift can be taken as approx- 
imately constant, A les (E) ~ A ICS (E = 0). The energy 
width T(E) = 2Tr\( ( j )b \HQp\* p -(E))\ 2 is proportional to 
k in the limit for k J. 0, which is a consequence of the 
surviving s-wave part of the scattering wave function. 
Therefore it is usually assumed that T(E) ~ 2Ck, with 
C a constant that characterizes the coupling strength be- 
tween P and Q |3J|. The resulting expression for the S 
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matrix 

S(k) = S P (k) ( 1 n 2lCk J (16) 

will be referred to as the single-resonance approximation. 

The direct part of the S matrix is related to the direct 
part of the scattering length as Sp(k) — exp[— 2ifcabg], 
and the second term describes the resonant behavior due 
to coupling to Q-space. In the limit k [ this single- 
resonance expression for the S matrix gives the dispersive 
behavior of Eq. 

However, the single-resonance approximation has a 
limited range of validity. More specifically, it is based 
on the assumption that the scattering wave function 
l^/p) can be approximated by its low-energy s-wave limit, 
which is proportional to k 1 ^ 2 (also known as Wigner's 
threshold law). If the P-channel has a low-energy res- 
onance, such as a (nearly-)bound state (with e^ s = 
— fi. 2 K 2 /2/i) close to threshold, this approximation breaks 
down when \E/e p cs \ <C 1 is not satisfied. 

In the following we show that the assumptions about 
the trivial energy dependence of Sp(k) and A(E) are not 
valid in the case the P-channel is nearly resonant. We 
derive how the full energy dependence can be taken into 
account analytically. In Sec. IIVI we discuss the Mittag- 
Leffler expansion of the resolvent operator \E + — Hpp]^ 1 
based on Gamow resonance states, which takes these P- 
channel resonances into account. In Sec. we derive an- 
alytical expressions for the energy shift and width based 
on this expansion. 



IV. OPEN-CHANNEL RESONANCES 

Without coupling to the closed-channel subspace Q, 
the scattering properties in the open-channel subspace V 
are governed by the uncoupled Schrodinger equation [E— 
Hpp)\$>p) = 0, which has scattering solutions |\&p). We 
assume that the interatomic potential is of finite range, 
i.e., it can be neglected beyond some radius Rq and is of 
the form 

t t ,x ( V PP {r) for0<r<i? , , . 

Vpp W = { for r > R . (17) 

The potentials used in the description of weakly inter- 
acting cold alkali gases have the asymptotic form V(r) ~ 
—Ce/r 6 . If we take Rq such that V(Rq) is small com- 
pared to the asymptotic kinetic energy E, the finite-range 
approximation can be applied retaining high accuracy 

The direct part of the S matrix is determined by 
the uncoupled P-space scattering solutions. In the low- 
energy limit the scattering properties for single-channel 
problems are determined by the s-wave part of the scat- 
tering wave function ipp (r) = (r\ipp). If there are no low- 
energy potential resonances, the scattering length ab g of 



the uncoupled Hpp problem encapsulates the direct part 
of the scattering process in the low-energy domain. 

However, if Hpp possesses a low-energy potential res- 
onance this is not true in general. The Sp matrix 
contains a resonant feature and becomes strongly en- 
ergy dependent near the collision threshold. The di- 
rect part of the S matrix takes the form Sp = S P S S P S , 
where the background part of the S matrix, described 
by S p — exp[—2ika p ], summarizes the effect of the P- 
channcl potential on the scattering solutions, excluding 
the resonant part. The resonant part is described by 
S P B . Note that we added a superscript P to distinguish 
the P-channel resonance from the Feshbach resonance, 
which is induced by coupling to the Q-channel. 

The resonances and bound states of the P-channel in- 
teraction potential correspond to the poles of Sp. More 
specifically, the energy dependence of the resonant part 
S p cs is determined by the anal ytic al properties of the (P- 
channel) Jost function !F{k) |4J|. The Jost function is 
related to Sp matrix as Sp[k) = fe)/J 7 (fe), and the 
zeros of the Jost function correspond to the poles of the 
Sp matrix. 

The Jost function T(k) is an entire function in k for 
finite-range potentials. In the general case, the zeros of 
T(k) in the upper halfplane are only located at the imag- 
inary k axis, are finite in number, and correspond to the 
bound states in V space. Note that since E < in this 
case, the V subspace is in fact also energetically closed. 
There is a countable, infinite number of zeros of !F(k) 
in the lower halfplane. The poles located in the third 
and fourth quadrant correspond to the resonances associ- 
ated with a non-zero angular momentum of the colliding 
particles, and are located symmetrically with respect to 
the imaginary k axis. The poles located at the negative 
imaginary axis are finite in number and correspond to 
the s-wave virtual states. Labelling the zeros by k n , the 
Jost function can be written as the product |44| 

°° / k \ 

n— 1 ^ n ^ 

where h(k) is a smooth function related to the range 
of the potential Tq. This expression is normalized such 
that it reproduces the correct behavior near threshold, 
F{k) — > 1 for k — > 0. Note that the physical parameter 
ro should not be confused with the cut-off radius Po, and 
typically r <C P - 

From the relation E = h 2 k 2 /2[i it follows that the 
mapping from k to E is two-to-one. As a result the single- 
valued functions of k, such as the Jost function and the 
S matrix, will be double- valued functions of the complex 
energy E. However, we can map the upper and lower half 
of the complex k plane to two different complex energy 
planes. The upper half (Im(fc) > 0) is mapped onto 
the so-called physical (or first) sheet, whereas the lower 
half (Im(fc) < 0) is mapped onto the non-physical (or 
second) sheet. If we consider the E plane as a two-sheeted 
Riemann surface |45| . the afore mentioned functions are 
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single- valued functions of the energy E. Both sheets have 
a branch-cut discontinuity on the real axis running from 
to oo. The two Riemann sheets are connected by the 
positive real .E-axis, where Im(fc) = 0. 

The bound states and resonances of a scattering sys- 
tem correspond to the poles of S(k). If we consider the 
scattering matrix as a function of the collision energy, 
which we will loosely write as S(E), there is a similar 
relation between the poles of S(E) on the two Riemann 
sheets and the bound and resonance states of a scatter- 
ing system. The bound states correspond to the poles of 
S(E) at the negative real E axis on the physical (or first) 
Riemann sheet. The I ^ resonances correspond to the 
poles of S(E) in the first and fourth quadrant of the com- 
plex E plane on the non-physical (or second) Riemann 
sheet. The s-wave virtual states correspond to the poles 
of S(E) at the negative real E axis on the non-physical 
Riemann sheet. 



Virtual states 

If the potential V effectively gets less attractive, the 
bound states of the potential will move towards the col- 
lision threshold. At some point the highest bound state 
will cross the collision threshold, and turn into a virtual 
state in the case of s-wave collisions. The correspond- 
ing transition of a bound-state pole into a virtual-state 
pole is schematically illustrated in Fig. |SJ A virtual state 
behaves much like a bound state in the inner region of 
the potential, but it is not a proper bound state as it be- 
haves asymptotically as e + ' K ' r . The important effect of a 
low-energy virtual state is that it introduces a resonance 
feature in the Sp matrix. 

In Sec. IVll we calculate the position of the poles of the 
Sp matrix in the case of 85 Rb, and find a virtual state at 
£virtuai/fc_B = —6.45 fj,K. This is close to threshold, since 
the typical energies in current experimental setups are 
of order microKelvin. We discuss how this virtual state 
can be accounted for in the following. The zero of the 
Jost function associated with the virtual state is labelled 
by ~iK vs , where n vs is a positive constant. Furthermore, 
all other poles of the Sp matrix are located far from the 
origin k = 0, and the contribution of these distant poles 
to the Sp matrix is summarized by a smooth function 
g(k). The Jost function can then be written as 



Im(<r) 



Im(£) 



T{h) = i + — g(k). 



(19) 



For real k it follows that the direct part of the S matrix 
is given by 



S P (k) = e~ Mka ^- 



2ikaf * K vs k 



k' 



(20) 



where the background factor exp[— 2ika bg ] summarizes 
the effect of all the non-resonant poles of the Sp matrix. 
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FIG. 6: Schematic illustration of the movement of the poles of 
the S matrix in the complex k plane (a) and the two-sheeted 
complex E plane (b). The solid E axes correspond to the 
first (or physical) Riemann sheet, while the dashed E axes 
correspond to the second (or non- physical) Riemann sheet. 
Both sheets have a branch-cut discontinuity on the real axis 
running from to oo. The two Riemann sheets are connected 
by the positive real axis, where Im(fc) = 0. Initially, the pole 
is located at the positive imaginary k axis, which corresponds 
to a bound-state pole on the negative real E axis on the first 
Riemann sheet. If the interaction potential effectively gets 
less attractive the pole moves towards k = 0, at some point 
crosses the origin, and turns into a virtual-state pole at the 
negative imaginary fc-axis. This corresponds to a virtual-state 
pole at the negative real _B-axis on the second Riemann sheet. 



The scattering phase 6(k) is related to the Sp matrix 
as Sp(k) = exp[2id(k)], and is evaluated as 



8(k) — —ka bg + arctan 



(21) 



The background part is related to the phase of g(k), and 
is linear in k in the low-energy limit .46]. The resonant 
part is related to the pole at k = —ik V s, and causes a 
'bump' in the scattering phase at low energies. Moreover, 
if the virtual-state pole gets closer to threshold (k vs — > 0), 
the scattering length ab g = a^ g — 1/ k V s will become more 
and more negative. 

To show that this virtual state has to be taken into 
account explicitly, Fig. [7| shows the scattering phase for 
85 Rb in the |/,m/) = |2, — 2} P-channel, without cou- 
pling to the Q-channels. The black dots represent the 
numerical results, which are obtained by solving the 
Schrodinger equation using the proper physical and state- 
of-the-art rubidium potentials. The solid line is ob- 
tained from the virtual-state expression, and the dashed 
line is obtained from the usual contact potential ap- 
proximation 8(k) — — arctan[fccibg] |35|. Let us stress 
again that Fig. shows the scattering phase for the P- 
channcl only, and the coupling to the Q-channel has 
been excluded from this calculation. For this particu- 
lar channel ab g = — 443 ao 14]. The virtual state con- 
tributes to the scattering length as — 1/k vs = — 562 ao, 
and a^ g = +119 ao is now of the order of the potential 
range ro- More details about the numerical calculation 
can be found in Sec. IVII 
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E(pK) 

FIG. 7: Scattering phase &(k) for 85 Rb in the \f,m f ) = |2, -2) 
spin channel. The black dots represent the numerical re- 
sults. The solid line is obtained from Eq. (1211 . The dashed 
line is obtained from the contact potential approximation 
S(k) = — arctan[fc<2b g ]. The scattering length of 85 Rb in 
this channel is db g = — 443 ao, in agreement with coupled- 
channels calculations. The background scattering length is 
a^g = +119 ao, and the virtual state is located at k VB — 
-k V8 = —1.78 ■ 10 _3 i [1/ao], which corresponds to the term 
— 1/kvs = — 562 ao in the expression for ab g . 



Comparing the virtual-state expression for S(k) with 
the numerical results, the agreement is excellent. If we 
compare the numerical results with the contact potential 
expression, it is immediately seen that this expression al- 
ready starts to deviate significantly at a few microKelvin. 
This indicates that the scattering length parameter only 
does not fully encapsulate the energy dependence of the 
scattering physics, and the P-channel resonance should 
be taken into account explicitly. 



Mittag-Leffler series 

The S matrix, T matrix, and the resolvent (or Green's) 
operator Gpp(E) = \E-Hpp\~ 1 have their poles in com- 
mon ^tJ- This suggests it is possible to expand the re- 
solvent in a Mittag-Leffler series , where the resolvent 
is written as a sum over the different pole contributions. 
The corresponding residues are related to the bound and 
resonance states, which are eigensolutions of Hpp. The 
scattering and transition matrix can then be written as 
a sum over the pole contributions as well. 

Eigensolutions associated with poles of the scattering 
matrix were first introduced by Gamow in his theoretical 
description of a-decay |48j . In the case of s-wave scatter- 
ing, the Gamow function fi n (r) = (r|J7„) associated with 
the pole k n of the scattering matrix, is defined as the 
solution to the radial Schrodinger equation that satisfies 



the boundary conditions 



0. 



^Q„(r)| r =R = ik n fl n (r)\ r=1 



(22) 



The radius R has to be chosen such that R > R , i.e., 
it is applied in the asymptotically free (V = 0) region of 
the interaction potential. 

The Gamow states behave asymptotically as f2 n (r) oc 
exp[ifc n r]. For the bound state poles k n = iK n , with n n a 
positive constant, the Gamow states are just the (prop- 
erly normalized) bound state wave functions. However, 
for the poles with Im(fc„) < the Gamow states expo- 
nentially diverge. Due to this exponential divergence, 
the Gamow states do not form an orthonormal basis for 
the V subspace. Defining a dual set of Gamow states 
as il® = fi*, the Gamow states do form a biorthogo- 
nal set in the sense that {fl^\n n >) — S nn >. Here S nn > is 
the Kronecker delta, and the inner product is defined by 
means of analytic continuation in k of the proper bound 
state eigensolutions to the resonance poles in the lower 
half of the complex /c-plane 0, • The normalization 
condition takes the form 



Jo 



2k n 



Q 2 n (R) = I. 



(23) 



Note that, as before, R has to be chosen such that R > 
Rq, in which case the normalization condition does not 
actually depend on the precise choice of R. 

The Gamow state \£l n ) is an eigenstate of Hpp and has 
an eigenvalue E n . The dual state |f2^) is an eigenstate 
of Hpp with the eigenvalue -E*. This can be written in 
Dirac notation as 



Hpp\n n ) 
(n°\H PP 



E n \£l n ), 



(24) 



From these last expressions it is immediately seen that 
the Hamiltonian Hpp is diagonal with respect to the 
biorthogonal set of Gamow functions. The matrix el- 
ements of Hpp can be evaluated as (Q.®\Hpp\Q, n i) — 

En' &nn' • 

Using the Gamow resonance states, the resolvent op- 
erator [E — Hpp]^ 1 can be expanded as a sum over 
eigenstates associated with poles of the scattering ma- 
trix. This type of expansion is known as the Mittag- 
Leffler expansion. The properties of the Mittag-Leffler 
expansion in terms of the Gamow resonance states 
have been extensively studied in the literature, see e.g. 

Ref. mm mm. 

We will therefore not give a deriva- 
tion of the Mittag-Leffler expansion, but refer the reader 
to the literature. We here just give the Mittag-Leffler 
expansion of the resolvent operator: 



i ~ |n n )(n?| 

E-Hpp ^2k n {k-k n )' 



(25) 



where n runs over all poles of the Sp matrix. For 
notational convenience, we have used units such that 



10 



fi = 2(1 = 1 so that k corresponds to the energy as 
E = k 2 . Note that E can be any complex energy in 
Eq. 125|) in principle. In Sec. we discuss how the 
Mittag-Leffler expansion can be used to find the energy 
dependence of A(E) analytically. 



V. INTERPLAY BETWEEN FESHBACH 
RESONANCES AND POTENTIAL 
RESONANCES 

Inserting the Mittag-Leffler series Eq. (|25[1 in the com- 
plex energy shift Eq. l|15fl results in 



A{E) = g (0 fc |g Q p|»n)R?|gpQ|^) 



n=l 



(26) 



The complex energy shift of the unperturbed bound-state 
in the Q-channel is thus expanded over the contribution 
of the various poles of the P-channel scattering matrix. 
The exact expression for the complex energy shift de- 
rived here is the most important difference with previous 
work on the interplay between single-channel and multi- 
channel resonances. 

We now illustrate the impact of a virtual state close to 
the collision threshold on the coupling between the spin 
channels P and Q. The full S matrix has the form S = 
S^ g S^ cs S^, s , and the poles of this S matrix correspond 
to the energies of the dressed (quasi-)bound states. In 
Sec. IIVI we already mentioned that in the case of 85 Rb 
there is a virtual state close to threshold. We assume that 
this pole is dominant, and the effect of all other poles of 
the P-channel can be summarized by a background term. 
The energy shift due to coupling then takes the form 



A(E) 



(27) 



The basis |fi n ) is energy independent, as is the Q- 
channel bound state \(j>b)- The coupling Hpq = Hq P 
is energy independent as well. It then follows that 
the energy dependence of A(E) is fully determined by 
the denominator in Eq. (|27J) . which is proportional to 
l/(fc-fc vs ). 

As the numerator of Eq. I|27[) is a real constant and 
fc vs = —in vs , we can write the complex energy shift as 



A(E) = 



-iA v 



(28) 



where A vs = -(4>b\HQp\Qr VS )(Q,^\HpQ\<t>b) is a positive 
constant. The extra minus sign in this constant shows 
up for virtual states due to the normalization condition 
of the Gamow states. 

For real and positive k (or equivalently, for real and 
positive E on the physical Riemann sheet) we can multi- 
ply the numerator and denominator with (k — iK vs ) and 



take the real and imaginary part of A(E): 



A re8 (£) 
T{E) 



2 ^vs 



k 2 



A VF ,k 



^vs 

(k 2 + K 2 VS ) 



(29) 
(30) 



For k = in on the positive imaginary axis (or equiva- 
lently, for real and negative E on the physical Riemann 
sheet) the energy shift is real- valued and given as 



A res (S) 



I A 



K vs (k + K V s) 



(31) 



The energies of the dressed states are given by the poles 
of the full S matrix, and can be found by solving 



(k + i Kvs ){E-e b (B) - A(E)) =0. 



(32) 



We see that in the presence of a virtual state pole close to 
threshold, the denominator of A YOS (E) gets close to zero 
for E — > 0. Therefore, the energy of the dressed molecu- 
lar Feshbach state ebind = £6 + A ros (P), strongly depends 
on the energy E close to threshold. For positive ener- 
gies, the energy shift A res (P) and width T(E), related to 
the quasi-bound Feshbach state, depend strongly on the 
energy as well. 

In order to give an accurate description of the dressed 
(quasi-)bound Feshbach state, it turns out that in the 
case of 85 Rb the influence of the highest P-channel bound 
state on the energy shift and width should be taken into 
account as well (see Sec. IVIfl . As can be seen in Fig. ITU1 
this P-channel bound state and the Q-channel bound 
state have a broad avoided crossing. In the inset, the ef- 
fect of the virtual state close to threshold is visible. How- 
ever, on a larger scale it is clear that the (quasi-)bound 
state has still not converged to the bare Q-state energy. 

The expressions for the energy shift and width includ- 
ing both pole contributions are given in the following. 
The pole associated with the highest P-channel bound 
state is denoted by fcbs = *^bs- For real and positive k 
(or equivalently, for real and positive E on the physical 
Riemann sheet) we take the real and imaginary part of 
A(E): 



T(E) 



2 ^vs 



A V vk 



l4u 



A hs k 



(33) 
(34) 



Here A\, s is a positive constant related to the coupling 
of the P-channel bound state to the Q-channel bound 
state. For k = in on the positive imaginary axis (or 
equivalently, for real and negative E on the physical Rie- 
mann sheet) the energy shift operator is real-valued and 
given as 



A res (£) 



2 Ay 



\A 



bs 



^vs (K + Kbs) 



(35) 
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The constants A vs and Ab s can be found from input of 
coupled-channels calculations, or equivalently from mea- 
surement (see Sec. IVIjl . Note that the contribution of 
the virtual state to the energy shift is negative valued, 
and the contribution of the proper bound state is posi- 
tive valued. Let us stress that there are no free param- 
eters in this model. The model is fully characterized by 
only a few parameters, which are determined by physical 
quantities directly related to the true interaction poten- 
tials. These parameters can be extracted from numeri- 
cal coupled-channels calculations, or directly from exper- 
imental measurements. More details are given in Sec. lVIl 



VI. NUMERICAL METHOD AND RESULTS 

The numerical results in this section are based on 
coupled-channels calculations [54| for rubidium, based 
on the most recent knowledge of the interaction poten- 
tials 0, |2fj . These calculations take the coupling be- 
tween the relevant spin channels into account, i.e. P and 
Q are coupled. In order to find the physical properties of 
the open channel only, we 'turn off' the coupling between 
P and Q. This allows us to calculate the single-channel 
Sp matrix, and find its poles and the background scatter- 
ing length. In order to find the other model-parameters, 
we perform a calculation where the coupling between P 
and Q is again taken into account. In principle, only two 
data points are needed to determine the constants A vs 
and A^ s . 

In Fig. |S1 the Sp matrix is shown for negative energies 
on the non-physical Riemann sheet, i.e. Im(fc) < 0, for 
85 Rb in the 2, —2) hyperfine channel. The pole of the Sp 
matrix on this sheet is located at e V irtuai/fc_B = —6.45 ^K. 
This indicates that the P-channel interaction potential 
has a virtual state at this energy. The phase of the Sp 
matrix for positive energies (on the physical sheet) was 
already shown in Fig. 

In order to compare the results of our model with the 
numerical coupled-channels calculation, we perform a full 
calculation over a large energy and magnetic field range 
around threshold and resonance. For positive energies 
with respect to the P-channel threshold, the scattering 
phase 6(E, B) is calculated. At each energy value the 
derivative of S(E, B) with respect to B is calculated. Ac- 
cording to Eq. (|14fl the derivative is given by the following 
expression: 

dS(E,B) A^T(E)/2 

dB [ e Q(B) + A rcs (P) - EY + [T{E)/2f ' 

Here A/x mag is the relative magnetic moment of the bare 
Q-channel bound state with respect to the P-channel 
threshold. This derivative function is the well-known 
Lorentz curve. The center of the Lorentz curve is given 
by the condition E = e£ (B) + A res (P), which determines 
the location of the dressed quasi-bound Feshbach state. 
The energy width of this state equals the width of the 




FIG. 8: Sp matrix for negative energies on the non-physical 
Riemann sheet for 85 Rb in the |2, —2) hyperfine channel. The 
pole of the Sp matrix on this sheet is located at e v irtuai/fcs = 
-6.45 (iK. 



Lorentz curve. We can thus calculate the position and 
width of the quasi-bound Feshbach state as a function of 
the collision energy. 

For negative energies with respect to the P-channel 
threshold, the coupled-channels Schrodinger equation is 
integrated outward, starting in the inner region of the 
interaction potential to some matching radius r m . The 
Schrodinger equation is also integrated inward from some 
r max to the matching radius r m . The boundary condi- 
tions at r max are such that the solution asymptotically 
vanishes and is physically acceptable. If the two solu- 
tions obtained in this way are linearly dependent, the 
corresponding energy is an eigenenergy. At this energy, 
the full coupled potential has a bound state. This calcu- 
lation is repeated as a function of magnetic field, and thus 
gives the field-dependent position of the dressed molecu- 
lar Feshbach state. 

The energy width of the dressed quasi-bound Feshbach 
state as calculated with a coupled-channels method is 
shown in Fig. (black dots). The dotted line shows the 
energy width according to our model, where only the P- 
channcl virtual state is taken into account. The constant 
A vs is determined by comparing Eq. (|3L)|) with a single 
low-energy data point. We see that for low energies the 
virtual-state expression and the coupled-channels data 
agree very well. However, for energies larger than roughly 
10 microKelvin, the virtual-state expression starts to de- 
viate. 

The solid line shows the energy width according to our 
model, taking the highest bound state of the P-channel 
into account as well. The energy width is now accu- 
rately described by our model for a large energy domain. 
We obtain the two parameters A vs and Ab s from a fit of 
Eq. I|34|) to two data points. 
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FIG. 9: Energy width of the dressed quasi-bound Feshbach 
state for 85 Rb in the |2, —2) hyperfine channel. The dot- 
ted line shows the energy width according to Eq. B3UP . where 
only the virtual state, with k vs = 2.54 ■ 10 -3 K 1 ^ 2 is taken 
into account, and A VB = 1.94 ■ 10~ 8 K 2 . The solid line shows 
the energy width according to Eq. (1341 . where the high- 
est bound state of the P-channel is included as well, with 
ft bs = 0.103K 1/2 . In this case ^ V8 = 1.92 • 10" 8 K 2 and 
,4b s = 1-26 • 1(T 5 K 2 . 



We now insert these parameters into the expressions 
for the energy shift A res (P) to describe the energy of 
the dressed (quasi-)bound Feshbach state. The result 
is shown in Fig. ^3 where the black dots indicate the 
coupled-channels results and the solid line is obtained 
from our model. In Fig. 1111 we zoom in closer to the 
P-channel threshold. 

The unperturbed (or bare) bound state in the Q- 
channel subspace is dressed by the coupling to the P- 
channel subspace. This induces an avoided crossing with 
the highest P-channel bound state. The avoided cross- 
ing is broad in the sense that, even though the un- 
perturbed P-channel bound state is located at roughly 
ebound/fcs ~ — 10 mK, close to the P-channel threshold 
the dressed state still has not converged to the bare Q- 
channel bound state. 

The P-channel virtual state is not located at the phys- 
ical energy sheet, and there is no avoided crossing of the 
usual kind between the dressed Q-channel (quasi-)bound 
state and the virtual state. However, the virtual state 
is located close to the collision threshold and induces a 
strong threshold effect. This threshold effect dominates 
the behavior of the molecular binding energy near the 
collision threshold, and has to be taken into account ex- 
plicitly. In our model we take the relevant P-channel 
bound and/or virtual states into account analytically. 
From these figures it is immediately seen that our model 
agrees perfectly with full coupled-channels calculations 
for a very large energy domain. The binding energy of 
the dressed molecular state that has been measured in 
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FIG. 10: (Color online) Energy of the dressed (quasi-)bound 
Feshbach state. The black dots indicate the coupled-channels 
data. The thick solid line indicates the energy according 
to our model. The thin solid line is the energy of the un- 
perturbed Q-channel bound state, which around threshold 
is given by the linear expression ef(B) = A//" ag (_B — Bo), 
with Afi mae = -1-75 ■ 10" 4 K/G = -3.64MHz/G and B = 
160.1 G. 
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FIG. 11: (Color online) Same as in Fig. 1101 but now for en- 
ergies closer to the P-channel threshold. 



Ref. 0, Q] is described analytically with high precision. 
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VII. OTHER APPROACHES 

In this section, we compare our model with some other 
approaches commonly used in the description of Fesh- 
bach resonances. A model which is conveniently used 
in many-body theories is the contact potential (or zero 
range potential). In this approach the real interaction 
potentials are replaced by deltafunctions or pseudopo- 
tentials, proportional to the s-wave scattering length a. 
In the vicinity of a Feshbach resonance the dispersive 
formula Eq. (JJ is used. The scattering matrix takes the 
form El 



S(k) 



1 - ika(B) 
l + ika(B)' 



(37) 



The molecular binding energy is determined by the pole 
of S(k), and is given as et,md = —h 2 /(2fia 2 (B)). In 
Fig. 1121 we compare the resulting binding energy (dot- 
ted line) with the coupled-channels results (black dots). 
It is clear that close to the resonance magnetic field Bo, 
the contact model binding energy agrees quite well with 
the coupled-channels binding energy. Further away from 
the resonance, the contact model energy starts to deviate 
significantly from the exact binding energy. 

In Eq. Ijl6|l the S matrix for the single-resonance Fesh- 
bach model is given. The pole of this S matrix gives the 
energies of the dressed (quasi-)bound state. This energy 
is indicated by the dashed-dotted line in Fig. ^] It can 
be clearly seen, as expected, that the single-resonance 
model fails already relatively close to threshold. 

Another commonly used model originates from the ef- 
fective range approach. In the contact model the scatter- 
ing phase S(k) is approximated as S(k) = — arctanffca], 
but this approximation already breaks down at low en- 
ergies. In the effective range model a second term is 
included to describe 5(k) further away from the collision 
threshold. The scattering matrix is given as 



S(k) = 



-l/a(B) +r a (B)k 2 /2 + ik 
-1/ 'a(B) + r (B)k 2 /2 - ik' 



(38) 



where r${B) is the effective range parameter 44]. The 
molecular binding energy can again be found by solving 



^- ) +lr (B)k^-ik = 0, 



(39) 



and is shown in Fig. ^] (dashed line). We determined 
ro(-B) and a(B) simultaneously from a coupled-channels 
calculation, by fitting to the scattering phase for two dif- 
ferent energy values, for several values of the magnetic 
field. Although the field-dependent effective range model 
agrees reasonably well with the coupled-channels binding 
energy for magnetic fields close to Bo, it breaks down for 
fields larger than B « 164 Gauss and cannot be applied 
further away from threshold. 

Summarizing, the contact model, single-resonance 
Feshbach model, and effective range model are low- 
energy approximations of the exact scattering matrix. 



These approaches give reasonable agreement with the 
coupled-channels molecular binding energy close to the 
resonant magnetic field Bo, but for energies further away 
from threshold these descriptions give poor agreement. 

Our model has several important advantages compared 
to the models discussed here. First of all, our model is 
directly based on the underlying physics of the interplay 
between potential and Feshbach resonances. The con- 
tact potential model is based on the assumption that the 
scattering phase can be replaced by its low-energy limit, 
8{k) — — arctan[faz(-B)]. As we have seen, the poten- 
tial resonance and Feshbach resonance introduce addi- 
tional energy-dependence of the scattering phase shift. 
The single-resonance Feshbach model would give an ad- 
equate description for the situation where ab g is on the 
order of the range of the potential. However, when o,b g is 
large compared to this range, this model breaks down as 
well. The effective range approach gives a better descrip- 
tion of the binding energy compared to the two previous 
models. However, it still breaks down at some point, 
and moreover, it does not give much physical insight in 
the mechanism behind the additional energy dependence 
of various cold-collision properties. Another clear disad- 
vantage of the effective range approach is that we have to 
calculate ro (B) numerically as a function of the magnetic 
field using a coupled-channels method. In our model the 
binding energy is given by a simple analytical formula, 
which gives an excellent description even for magnetic 
fields far from Bo and energies of order milliKelvin. Our 
model is fully characterized by only a few parameters, 
which can be extracted from coupled-channels calcula- 
tions, or directly from measurements. 



Two resonances in Q-space 

Another approach that has been proposed in Ref. [5j| 
is to use a double-resonance parametrization of the scat- 
tering matrix within the Feshbach projection formalism. 
The second resonance introduced should account for the 
influence of the potential resonance on the properties of 
the Feshbach resonance. 

The Feshbach projection formalism can be used to 
describe potential resonances, albeit in a rather indi- 
rect manner. In a paper by Domcke |33j]. the Feshbach 
projection-operator approach is used to describe poten- 
tial resonances in scattering systems. The projector onto 
the resonance states is defined as Q = Y^n=i\^ > n)(4'n\> 
where the set of functions is an arbitrary or- 

thonormal set of square- integrable functions. The formal 
requirement on the projection operator Q is given as 



(r\Q\1>) 







for any \ip). 



(40) 



The formalism leads to a decomposition of the scattering 
matrix into a resonant and a non-resonant part, where 
the resonances are all contained in Q space if N is cho- 
sen large enough. The choice of the states \<j> n ) to be 
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introduced above, and the total scattering matrix takes 
the form loal 



-100 



FIG. 12: (Color online) Comparison of coupled-channels cal- 
culations of the binding energy (black dots) with the con- 
tact model (dotted line), the single-resonance Feshbach model 
(dashed-dotted line), the effective range model (dashed line), 
and our virtual state model (solid line). Although for the 
magnetic fields shown the effective range model and our model 
give comparable results, our model gives an analytical expres- 
sion for the binding energy. This is not the case for the effec- 
tive range model, where a field-dependent parameter ro(B) 
has to be taken into account that has to be calculated numer- 
ically. 



used is arbitrary, however, and the correspondence be- 
tween these states and the real resonance states (such as 
the Gamow states) is not clear. More specifically, the 
Gamow states in general do not satisfy relation Eq. H4U|) 
and cannot be used to construct an operator onto the 
real resonance states within this formalism. 

However, as Domcke argues, increasing the number of 
states used in the projector Q by one, seems to remove 
exactly one resonance pole from the background part of 
the scattering matrix. One could therefore hope that 
adding an appropriate bound state \4>i), the virtual state 
resonance due to the Gamow resonance state |n v i r t ua i) 
can be indirectly included in the Q subspace. This is the 
approach followed in Ref. [55| , in order to account for the 
interplay between a Feshbach resonance and a potential 
resonance in 6 Li collisions. We will refer to this as the 
double-resonance approach. 

Although the bound state introduced to account for 
the second resonance is not clearly linked to the P- 
channel virtual state that gives rise to the potential reso- 
nance, the double-resonance approach is mathematically 
equivalent to our model under some constraints. We will 
show this in the following. 

Introducing two Q-channel bound states \<fii), one to 
account for the Feshbach resonance (i = b) and one to 
account for the virtual state resonance (i = v), the direct 
part of the scattering matrix does not contain a resonant 

-2ikaQ, 



feature anymore and takes the form S£ g = exp[ 
where is on the order of the range of the interaction 
potential. The Q subspace contains the two resonances 



S(k) — S^ g 



1 



2ik(C b A v + C v A b ) 
A b A v + ik{C b A v + C v A b ) 



(41) 



Here Aj = E — ef — A, is the detuning of the dressed state 
i with energy ef + Aj. The decay width of the dressed 
state i is given by = 2Cjfc. Note that the energy shift 
of both resonance states is approximately constant, and 
the decay width scales with k according to the Wigner 
threshold law. This is a direct consequence of the removal 
of the low-energy resonance from the direct channel by 
including the state \<p v ) into the Q subspace. 

The poles of the scattering matrix determine the (com- 
plex) energy of the dressed resonance states. In the 
double-resonance approach these poles are determined by 
the condition 



AfcA^ + ik(C b A v + C v A b ) = 0. 



(42) 



In our model the poles are determined by the condition 

(k + i Kvs )(E - a® - A(E)) =0, (43) 

if we only take the P-channel virtual state and the Q- 
channel bound state into account. In other words, we 
neglect the avoided crossing with the highest P-channel 
bound state. Under the constraint that the dressed 
bound state that removes the potential resonance from 
the direct channel is located far away from the collision 
threshold, i.e., A„ ~ — e v — A„, the two approaches are 
equivalent if the following relations are satisfied: 

A„ 



C v 

c b 



(44) 
(45) 



These relations show how the double-resonance parame- 
ters are related to the position of the virtual state (de- 
scribed by Kys), and the zero-energy shift of the dressed 
Feshbach resonance state, Ab = A res (P = 0). Note that 
the energy E always has to be negligible compared to 
e v + A„, otherwise the double-resonance S-matrix will 
introduce a non-physical energy-dependence in the scat- 
tering phase-shift. 

The double-resonance model does give an equivalent 
description of the scattering process, and can be used 
to parameterize the scattering matrix and/or scattering 
length. The link between the properties of the open- 
channel resonance and the second resonance introduced 
in the double-resonance model is not really clear. Our 
model has the important advantage that it is directly 
related to the underlying physics giving rise to the P- 
channel resonance. 



VIII. CONCLUSIONS 

In this paper, we have derived an analytical model that 
describes the cold-collision properties of two interact- 



15 



ing particles near a Feshbach resonance with large back- 
ground scattering length. The large background scatter- 
ing length results from an open-channel resonance near 
threshold, and this resonance has to be treated explicitly. 
The open-channel scattering is included in the Feshbach 
theory of resonances via a contribution from its poles of 
the S matrix, and its non-resonant open-channel back- 
ground scattering length. Here the latter corresponds to 
the true range of the potential. As an example, we study 
the B = 155 G Feshbach resonance of 85 Rb, and we 
show that our model compares excellent with numerical 
coupled-channels calculations in a large range of energies 
around threshold. 

Our model offers a simple physical picture for the 
understanding of threshold effects around resonance. 
It explains how the energy-shift between the dressed 
(quasi-)bound Feshbach state and the bare molecular 
state can be found. In the literature, several other models 
have been proposed to describe the formation off quasi- 
bound) molecules (see for instance 0,0,|Ii||53) We 
discuss in the following how our model compares to some 
other models. 

Several ideas have been proposed to describe the for- 
mation of (quasi-bound) molecules. In some descriptions 
(c.f. [HI,[l3,ISa]) the threshold and (bare) molecular state 
are regarded as a two-level system, with an effective cou- 
pling that accounts for the coupling between the open 
and closed channels. This results in a Landau-Zcncr 
crossing between the threshold state and the molecu- 
lar state. As can be seen from Fig. [SJ a Landau-Zener 
crossing does not correctly describe the dressed molec- 
ular state in case the open channel has a virtual state 
near threshold. The difference is striking, since the quasi- 
bound state in the continuum starts from threshold at Bq 
(where a is infinity), and then curves towards the bare 
Feshbach state. In the Landau-Zener model the quasi- 
bound state curves towards threshold at the other side 
of the bare Feshbach state. In case there is a real bound 
state in the open channel, the behavior of the dressed 
quasi-bound state appears to be qualitatively in agree- 
ment with the Landau-Zener crossing (c.f. Fig.QJ). How- 
ever, this is only the case since there is a real avoided 
crossing below threshold between the Feshbach state and 
the open-channel bound state. Threshold effects (visible 
in the inset) are not accounted for in the Landau-zener 
crossing. 

In other descriptions (c.f. |5j|) the energy levels of an 
interacting pair of atoms in a harmonic trap are calcu- 



lated. Here the lowest state behaves as the Feshbach 
bound state in the free-particle case for a small and posi- 
tive, and connects to the lowest harmonic oscillator state 
for a small and negative. All other trap states start on 
one side of the resonance (a > 0) in a particular trap 
state n, and are connected on the other side of the res- 
onance (a < 0) to a trap state n + 1. There is a strong 
field dependence of the eigenenergies for field values close 
to resonance, and the bare trap-levels can be regarded 
as having avoided crossings with the quasi-bound state. 
This can be understood from analytical models of two 
trapped atoms [|| EH, which calculate the energy 
shift of the trapped atoms due to interactions, using the 
scattering phase as a boundary condition for the trap- 
wavefunction of the particles. In this way, it can be seen 
that the level-energy is strongly field dependent when it 
coincides with the quasi-bound Feshbach state in the free- 
atom case. In the situation of a virtual state in the open 
channel (c.f. Fig|SJ), the quasi-bound state energy changes 
its slope in a region close to threshold. The avoided cross- 
ings between bare trap-levels and quasi-bound Feshbach 
state should therefore follow this behavior. 

The interplay between an open-channel resonance and 
a Feshbach resonance is usually not taken into ac- 
count in the description of resonance many-body sys- 
tems EH H il E3- The open-channel will 
have an important effect in a regime where the energies 
of the system are of order of or larger than the energy 
of the open-channel virtual state. A clear example where 
this effect is visible is the Ramsey-fringe experiment with 
85 Rb 0113, where the range of binding energies is so 
large that the bound state cannot be properly described 
by a single-resonance model. 
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